Formation of Multipartite Entanglement Using Random Quantum Gates 
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The formation of multipartite quantum entanglement by repeated operation of one and two qubit 
gates is examined. The resulting entanglement is evaluated using two measures: the average bi- 
partite entanglement and the Groverian measure. A comparison is made between two geometries 
of the quantum register: a one dimensional chain in which two-qubit gates apply only locally be- 
tween nearest neighbors and a non-local geometry in which such gates may apply between any pair 
of qubits. More specifically, we use a combination of random single qubit rotations and a fixed 
two-qubit gate such as the controUed-phase gate. It is found that in the non-local geometry the 
entanglement is generated at a higher rate. In both geometries, the Groverian measure converges 
to its asymptotic value more slowly than the average bipartite entanglement. These results are 
expected to have implications on different proposed geometries of future quantum computers with 
local and non-local interactions between the qubits. 

PACS numbers: 03.67.Lx, 89.70.-l-c 



I. INTRODUCTION 

Quantum algorithms such as Shor's factoring algo- 
rithm [l| and Grover's search algorithm 0, Q exemplify 
the potential speedup offered by quantum computers. Al- 
though the origin of this speed-up is not fully understood, 
there are indications that quantum entanglement plays a 
crucial role [1, [1, @1- Therefore, it is of interest to quan- 
tify the entanglement of the quantum states of multiple 
qubits that appear in quantum algorithms, using suitable 
measures. The special case of bipartite entanglement has 
been studied extensively in recent years and was estab- 
lished as a resource for quantum teleportation. The en- 
tanglement of pure bipartite states is completely deter- 
mined by their Schmidt coefficients. Thus, any measure 
of bipartite entanglement must be a function of these co- 
efficients. For example, the von Neumann entropy is ex- 
pressed in terms of the eigenvalues of the reduced density 
matrix, which are the squares of the Schmidt coefficients. 
For mixed bipartite states, several measures were pro- 
posed, namely the entanglement of formation and the en- 
tanglement of distillation 0, Q- In particular, for states 
of two qubits an exact formula for the entanglement of 
formation was obtained [1, ■ The more general case of 
multipartite entanglement, in a register of n > 2 qubits, 
is not as well understood, partly because no analogue 
of the Schmidt decomposition was found for multipartite 
systems. Axiomatic considerations have provided a set 
of p ro p erties that entanglement measures should satisfy 
pU . [T^ . [isl . [l^ . These properties include the requirement 
that any entanglement measure should vanish for product 
(or separable) states; it should be invariant under local 
unitary operations and should not increase as a result of 
any sequence of local operations complemented by only 
classical communication between the parties. Quanti- 
ties that satisfy these properties are called entanglement 
monotones. These properties provide useful guidelines 
in the search for entanglement measures for multipartite 
quantum states. Entanglement measures based on met- 
ric properties of the Hilbert space [ll|, [l^, [ill ^^'^ '^^ 



polynomial invariants |16l . |l7[ were proposed and shown 
to satisfy these requirements. Although some measures 
have been studied extensively, the connection between 
such measures and the efficiency of quantum algorithms 
remains unclear. 

The common framework of quantum computation is 
based on a universal set of one and two qubit gates 
[isl . [igj . Repeated operations of these gates enable to 
obtain any desired quantum state of the register. To 
generate multipartite entanglement encompassing all the 
qubits in a register, one applies a series of two qubit gates 
between pairs of qubits, in addition to single qubit gates. 
One can visualize this process as a network, in which the 
qubits are represented by nodes and two qubit operations 
are represented by edges that connect the corresponding 
nodes. Quantum entanglement is transitive in the sense 
that if we apply a gate that entangles qubits i and j and 
then a gate that entangles qubits j and k, this typically 
also gives rise to entanglement between qubits i and k. 
By entanglement between qubits i and k we mean that 
any two parties, each including one of these qubits, can- 
not be in a product state with each other. Thus, in or- 
der for the entanglement to encompass all the qubits in 
the register, the nodes associated with any pair of qubits 
must be connected, either directly or indirectly. 

In recent years, several possible implementations of 
quantum computers have been proposed. In some of 
these schemes the interactions between qubits are non- 
local in the sense that two-qubit gates may apply between 
any pair of qubits. Other schemes are based on a rigid 
geometry in which only nearest neighbor qubits may in- 
teract. In particular, they may be arranged in a one 
dimensional chain, as in solid state quantum computers 
using quantum dots [20| . In a one-dimensional geometry, 
two qubit gates may apply on only ti — I or n pairs of 
qubits (depending on the boundary conditions), out of 
the n{n— l)/2 possible pairs. It is interesting to examine 
to what extent the geometrical restriction reduces the ef- 
ficiency of the quantum computer. This limitation was 
addressed in the context of certain quantum algorithms 
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In this paper we consider the rate of formation of mul- 
tipartite entanglement in a register of n qubits by re- 
peated operations of one and two-qubit gates. Starting 
from a product state, we repeatedly apply a combination 
of random single qubit rotations and a fixed two-qubit 
gate on pairs of qubits. Two geometries are considered: 
a non-local geometry in which any pair of qubits may 
interact with each other and a one dimensional chain, in 
which two qubit gates may apply only between nearest 
neighbor qubits. The resulting entanglement is evaluated 
using two measures: the average bipartite entanglement 
[27l . \2§i . [29j and the Groverian measure [s^l ■ It is found 
that the non-local scheme is more effective, namely it 
requires fewer steps to produce a certain level of multi- 
partite entanglement compared to the local scheme. 

The paper is organized as follows. The entanglement 
measures used in this paper are briefly presented in Sec. 
II. The entanglement generating schemes, with local and 
non-local geometries are shown in Sec. III. The simula- 
tions and results are presented in Sec. IV, followed by a 
discussion in Sec. V and a summary in Sec. VI. 



II. ENTANGLEMENT MEASURES OF 
MULTIPLE QUBITS 

A. The Average Bipartite Entanglement Measure 

Consider a quantum register of n qubits in a pure state 
\ip). The bipartite entanglement between a given qubit 
and all the other qubits is given by a single parameter 
[Slj . namely, all the entanglement monotones for such 
partition are the same up to a monotonic function. One 
of these monotones is the largest eigenvalue P of the re- 
duced density matrix, p, of the single qubit. This eigen- 
value is also the square of the largest Schmidt coefficient. 
Another monotone is Tr((0^) = + (1 — P)^, which is 
a monotonically increasing function of P in the relevant 
interval 1/2<P<1. To evaluate the entanglement in 
a register of n qubits we average this measure over all 
n choices of the single qubit, where pi, i = 1,2, 
is the reduced density matrix of the ith qubit, obtained 
by taking a partial trace over all the other qubits. Af- 
ter a suitable normalization and shift, one obtains the 
entanglement measure [13, [11] 



g(^)^2--^Tr(p2). (1) 

i=l 

This measure was used in order to evaluate the entanglc- 
ment generated by repeated operations of random gates 
[1^1 Q(V') is an entanglement monotone only for the 
most refined partition, where each qubit is considered as 
a separate party. It is thus commonly considered as a 
measure of multipartite entanglement. 



In fact, the measure Q{ip) essentially quantifies bipar- 
tite entanglement, which is averaged over n different par- 
titions. In each partition, one party consists of a single 
qubit while the other party consists of n — 1 qubits. One 
may define other partitions of the register into two par- 
ties, where one party includes k qubits and the other 
party includes n — k qubits. For a given value of fc, the 
entanglement measure is obtained by averaging Tr(p^) 
over all possible partitions of this type. In a more gen- 
eral framework, the register can be partitioned into any 
number of parties between 2 and n [3, ■ 

B. The Groverian Measure 

Grover's algorithm performs a search for a marked ele- 
ment m in a search space D containing N elements [l, 0] . 
We assume, for convenience, that N = 2", where n is an 
integer. This way, the elements of D can be represented 
by an 7i-qubit register |a;) = \xi,X2, ■ ■ ■ ,Xn), with the 
computational basis states i = 0, ...,iV — 1. The 
meaning of marking the element m, is that there is a 
function f : D ^ {0, 1}, such that / = 1 for the marked 
element, and / = for the rest. To solve this search 
problem on a classical computer one needs to evaluate 
/ for each element, one by one, until the marked state 
is found. In the worst case, this requires N evaluations 
of /. On a quantum computer, where / can be evalu- 
ated coherently, Grover's algorithm, represented by the 
unitary operator Ug, can locate a marked element using 
only 0{^/N) coherent queries of / 0, [^. The algorithm 
is based on quantum oracle, with the ability to recog- 
nize the marked states 0, Starting with the equal 
superposition state, 

I'?) - E K)' (2) 

1=0 

and applying the operator Ug, one obtains the state 
Ug\v) ~ l*^) + 0{1/N), which is then measured. The 
success probability of the algorithm is almost unity. The 
adjoint equation takes the form (r/| — (m|J7G + 0{1/N). 
If an arbitrary pure state, \^), is used as the initial state 
instead of the state \ri), the success probability is reduced 
toP, = \{m\UGm^+0{l/N) OTP, = Mij)\^+Oil/N). 
The success probability is thus determined by the fidelity, 
between and \r]) fsol. [s^. 

Consider Grover's search algorithm, in which an arbi- 
trary pure state lip) is used as the initial state. Before 
applying the operator Ug, there is a pre-processing stage 
in which arbitrary local unitary operators Ui, U2, ■ ■ ■, 
Un are applied on the n qubits in the register. These 
operators are chosen such that the success probability of 
the algorithm will be maximized. The maximal success 
probability is thus given by 

PmaxW= max \{m\UG{Ui®...<»UnM)\\ (3) 
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which can be re-written as 

PmaxW = max |(</)|V')P, (4) 

where T is the space of all tensor product states of the 
form 10) = \(f>i) (g) . . . (g) \4in)- The Groverian measure is 
given by G{'4>) = \/l — Pmax (V') In the case of pure 
states, for which G{il)) is defined, it is closely related to 
the entanglement measure introduced in Refs. [ill. [l^. [l5| 
for both pure and mixed states and was shown to be an 
entanglement monotone. Based on these results, it was 
shown [sO] that: G(^/') > 0, with equality when |^) is a 
product state; G'(^) is invariant under local unitary op- 
erations and cannot be increased using local operations 
and classical communication. Therefore, G'('0) is an en- 
tanglement monotone for pure states. A related result 
was obtained in Ref. [3^, where it was shown that the 
evolution of the quantum state during the iteration of 
Grover's algorithm corresponds to the shortest path in 
the Hilbert space using a suitable metric. 

The Groverian measure was used in order to evaluate 
the entanglement generated by quantum algorithms such 
as Grover's algorithm and Shor's algorithm Js^]. It 
was also generalized to the case of mixed states [33| and 
to arbitrary partitions of the register [l^l . Here we use a 
different version of the Groverian measure, referred to as 
the logarithmic Groverian measure, given by 

G(V') = -ln(F,„axW). (5) 

This is an entanglement monotone, because the logarith- 
mic function is monotonically increasing, and it vanishes 
for product states, for which Pmax = 1- Unlike the mea- 
sure introduced in Ref. [l^l, which is restricted to the 
range of [0,1), the logarithmic Groverian measure may 
take values in the range < G{il}) < oo. This enables to 
better distinguish and compare between highly entangled 
states which involve a large number of qubits. 

Furthermore, the logarithmic Groverian measure ex- 
hibits the additivity property described below. Consider 
the state = |'!/)a)|V'b) where I^'a) and Ii/^b) are pure 
states of two different registers, of ua and ne qubits, re- 
spectively, and ua+ub = n. One can express the product 
state of n qubits in the form = |(/)a)|0b), where 
\4>a) and |0b) are product states of tt-a and ub qubits, 
respectively. Since is a tensor product of \iIja) and 
\4'b), we obtain that = (0a|'0a)(</'b|'0b)- Thus, us- 

ing Eq. (gl) we find that Pma.x{lp) = -Pmax(V'A)-Pmax(V'B)- 

As a result, the logarithmic Groverian measure satisfies 

G(^) = G(^a) + G(V'b). (6) 

The additivity enables to compare between the entan- 
glement in registers which include different numbers of 
qubits. 



III. THE ENTANGLEMENT GENERATING 
SCHEME 

To examine the rate of multipartite entanglement for- 
mation in quantum circuits we consider a scheme for 
the production of pseudo-random states of n qubits [1^ . 
The number of random amplitudes in a random state 
of n qubits is exponential in n. While the randomiza- 
tion of the state of each qubit alone requires only lin- 
ear resources, the formation of entanglement between 
them is more costly. Nevertheless, it can be achieved 
by a combination of random single qubit gates and a 
fixed two-qubit gate, which apply repeatedly on pairs 
of qubits. It was shown that this rotate- entangle-rotate 
scheme provides pseudo-random states with only polyno- 
mial resources (29l [39| . The random single qubit rotation 
is obtained using 

where at each step x, 71 and 72 are drawn from a 
uniform distribution over the unit interval [0,1). This 
parametrization provides operators drawn uniformly 
from the Haar measure For the two-qubit gate, we 
use the controUed-phase gate [2^ 

C/2 = 6^^'"-®'^% (8) 

expressed in the canonical decomposition form given in 
Ref. [3l|. This is a sensible choice because this operator 
can produce maximally entangled states. 

The initial state of the entanglement generating 
scheme is a random product state, {ipo), of n-qubits, 
where the state of each qubit is drawn from a uniform 
distribution on the Bloch sphere. Each iteration of this 
scheme consists of the following operations: 

1. Choose a random pair of qubits, i and j, from the n 
qubits in the register. In the non-local scheme, any 
pair of qubits may be chosen with equal probabili- 
ties. In the local scheme the qubits are arranged in 
a one dimensional chain with periodic boundaries 
and only nearest-neighbor pairs are chosen, with 
equal probabilites. 

2. Apply the controlled phase operator on qubits i and 
j- 

3. Apply random single qubit rotations, drawn uni- 
formly from the Haar measure [40| , on qubits i and 

j- 

The state obtained after t iterations, or time steps, is 
denoted by |'0t)- To examine the rate of entanglement 
formation and to compare between the two schemes, we 
evaluate the following three functions of Itpt)- (a) the 
fidelity F{'ipQ,'ipt) = KV'olV't)!; (t>) the average bipartite 
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entanglement Q{ipt)', and (c) the logarithmic Groverian 
measure G{ipt)- Unlike Q{tpt) and G{'ipt) which are en- 
tanglement measures, F{il'o,4't) is not an entanglement 
monotone. To put it on a common footing with the 
Groverian measure, we define 



(9) 



which is a monotonically decreasing function of the fi- 
delity and takes values between zero and infinity. 



IV. SIMULATIONS AND RESULTS 

Each simulation of the cntanglcmcnt-gcncrating 
scheme creates a series of quantum states \ipt), t = 
1,2,.... We examine the variation of K{ipQ,ipt) Qi'ipt) 
and G{ipt) vs. t. Due to the random nature of the op- 
erations, these three functions strongly fluctuate when 
evaluated for a single run of the entanglement generating 
scheme. To reduce the noise and elucidate the systematic 
trends, we perform a large number of runs using different 
initial states |V'o) and calculate the averages {K {ipo , tpt)) 
(Qi^t)) and (G(V't)> vs. t. 

In Fig. [T]we present {K{ipQ,ipt)) vs. t, averaged over 
2000 realizations, for the local scheme [Fig. [Ija)] and 
for the non-local scheme [Fig. [Ijb)]. The sharp increase 
of {K{ipo,tpt)) is explained by the fact that the fidelity 
is highly sensitive not only to two-qubit operations but 
also to single-qubit operations. Hence, the first few steps 
of the scheme can decrease the fidelity substantially. To 
explain the overshoot, notice that product states have 
smaller fidelity with each other than entangled states, 
since the fidelity of two product states is simply the mul- 
tiplication of the fidelities of the tensor components (all 
smaller then 1). In the beginning of the scheme, the 
register is still nearly separable so its fidelity with the 
(also separable) initial state is smaller than when it is 
completely entangled. After a sufficient number of steps, 
the state of the register is uncorrelated with the input 
state, and {K{ipQ^ ipt)) saturates. As the number of qubits 
in the register increases, it takes more steps to bring 
{K{ip{),'ipt)) to saturation. However, there is no signif- 
icant difference between the two geometries. 

The average bipartite entanglement, {Q{ipt)), is shown 
in Figs. [HJa) for the local scheme and in Fig. [D^b) 
for the non-local scheme. It increases more slowly than 
{K{ipQ, 4't))- The asymptotic value of {Q{tpt)) is the same 
for both geometries, but it is reached much faster in the 
non-local geometry than in the local, one dimensional 
geometry. The results for the one dimensionalgeome- 
try are consistent with those presented in Ref. [23[. The 
Groverian measure G'('0t) vs. t is shown in Fig. and 
Eljb), for the local and non-local geometries, respectively. 
Clearly, {G{'ipt)) converges to its asymptotic value more 
slowly than {Q{ip)). The asymptotic value of {G{ip)) is 
the same for both geometries but it is reached much faster 
in the non-local geometry. 
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FIG. 1: (Color online) The average of K(tpo,^t) = 
— In F{tpo, Ipt) vs. time in the local scheme (a) and in the 
non-local scheme (b) . The number of qubits in the register is 
n = 2, 3, . . . , 11 qubits (from bottom to top). 



The number of steps required for {K{i(jo, ipt)), {Q{ipt)), 
and {G{ipt)), to reach 90% of their saturation values, vs. 
n, is shown in Figs. IDJa), IHb) and ID^c), respectively. 
The results for {K{ipQ,ipt)) are almost identical in the 
local and non-local schemes. The time it takes {Q{ipt)), 
to reach 90% of its saturation values appears to be linear 
in n, for both geometries. However, the slope is lower 
for the non-local geometry, which means that the average 
bipartite entanglement builds up more quickly when non- 
local interactions are allowed. Similarly, in the non-local 
geometry the saturation time of {G{ip)) is linear in n. In 
the local geometry, the saturation time is longer and it 
deviates from linear dependence on n and is well fitted 
by a quadratic function of n. 

In the analysis above we focused on the averages of 
functions K, Q and G. Each data point was obtained 
by averaging over at least 2000 runs of the entanglement 
forming procedure described above. These averages were 
taken from the distributions f{K), f{Q) and /(G) for 
the quantum states \ipt) obtained after t iterations of the 
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FIG. 2: (Color online) The average of the entanglement mea- 
sure Qi^ilit) vs. t in the local scheme (a) and in the non-local 
scheme (b), for a register with 2, 3, . . . , 11 qubits, from bot- 
tom to top. Clearly, the entanglement builds up faster in the 
non-local scheme. 



procedure. As t increases, these distributions are found 
to approach the distribution obtained for random states 
of the register In Fig. 5 we present the probabihty 
densities /(G) of the Groverian measure G{ipt); for t ~ 
10 (dashed-dotted hne), 20 (dashed hne) and 100 (sohd 
hue) steps, for a register with n = 8 qubits. As the 
number of steps increases, these distributions approach 
the distribution /(G) obtained for pseudo-random states 
of the 8-qubit register (dotted Hne). 

These results indicate that building up bipartite as well 
as multipartite entanglement, using two-qubit gates is 
more efficient when non-local interactions are allowed. 
Furthermore, in both geometries it seems that multipar- 
tite entanglement, evaluated by G{'ip) builds up more 
slowly than the bipartite- like entanglement, quantified 
by QW. 




t 



FIG. 3: (Color online) The average of the Groverian measure 
G[4't) vs. t in the local scheme (a) and in the non-local scheme 
(b), for a register with 2,3, ...,11 qubits. The Groverian 
measure increases faster in the non-local scheme. 



V. DISCUSSION 

Quantum algorithms are typically designed under the 
assumption that two-qubit operations are possible be- 
tween any pair of qubits. However, some of the physical 
realizations of quantum computation allow only local in- 
teractions between neighboring qubits [l^]. In these re- 
alizations, non-local operations are achieved using s wap 
gates dll , which may involve a significant overhead [4 if . 
To put the results of this paper in the context of known 
physical realizations, we briefly review below several such 
realizations and specify whether they are based on local 
or non-local interactions. 

Ion-trap quantum computing devices (4^ . [4^ consist 
of laser-cooled ions, which are electromagnetically con- 
fined in ultra-high vacuum. The spin states of the ions 
are used as qubits. Single-qubit gates are preformed by 
selectively applying electromagnetic fields on the differ- 
ent ions. Two-qubit gates are realized through a global 
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FIG. 5: (Color online) The probability density /(G) of the 
Groverian measure G(V't) for the states obtained after f = 10 
(dashed-dotted line), 20 (dahsed line) and 100 (solid line) 
iterations of the entanglement forming procedure, for n = 8 
qubits. As t increases, the distribution converges to the one 
obtained for pseudo-random states (dotted line). 



FIG. 4: The number of steps required for the averages of 
K{iljo,tpt) (a), Q{ipt) (b) and G(V't) (c) to reach 90% of their 
saturation values in the local (+) and in the non-local (o) 
schemes vs. the number of qubits in the register. While the 
results for K{'il)o,ipt) are nearly identical in both schemes, the 
time required for Q{tpt) and G{tl)t) to converge is shorter in 
the non-local scheme. The solid lines are linear fits and the 
dashed line is a quadratic fit. 



phonon state, which makes use of the collective vibra- 
tional degree of freedom of the ions. The ground state 
and the first excited state of the phonon comprise a two- 
Icvcl system, which is used as an additional qubit. Us- 
ing carefully tuned laser-induced transitions, two-qubit 
gates can be applied between each ion and the global 
phonon. Specifically, a swap gate can be applied be- 
tween them, enabling to apply two-qubit gates between 
any pair of qubits. Ion traps thus provide a non-local 
implementation of quantum computation. Linear optics 
provides another non-local realization of quantum com- 
putation [H,!!!]. 

Quantum computation can also be performed using 
solid state devices [i^ HJ- Several proposals are com- 
prised of qubits based on superconducting Josephson 
junctions [4g|. One such proposal uses the two states of 
a superconducting single-charge box as a charge-qubit. 
The possibility of connecting all the qubits in parallel to 



a common LC-oscillator mode was suggested as a way 
to apply two-qubit gates [i^. This parallel connection, 
together with the ability to simultaneously control the 
Josephson coupling for each qubit, allows two-qubit in- 
teraction between any pair of qubits, giving rise to a 
non-local scheme. A different superconducting scheme is 
based on Josephson flux qubits, which interact through 
magnetic induction. Unlike the charge-qubit scheme, the 
flux-qubit scheme is a local one. Another implementa- 
tion is based on electron spins as qubits, confined to 
quantum dots [50| . Two-qubit gates arc applied using 
a time-dependent Heisenberg exchange coupling between 
adjacent dots, providing a local scheme for quantum com- 
putation. Solid state realizations may be most suitable 
for scalable implementation of quantum computations. 
Some of these realizations are based on local interac- 
tions. Local interactions appear form multiple-qubit en- 
tanglement at a lower rate than non-local interactions. 
However, the connection between the entanglement as 
quantified by the Groverian and related measures and 
the speedup offered by quantum algorithms is not yet 
clear. 



VI. SUMMARY 

We have studied the formation of multipartite quan- 
tum entanglement by repeated operation of one and two 
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qubit gates. The resulting entanglement was evaluated 
using the average bipartite and the Groverian measures. 
A comparison was made between two geometries of the 
quantum register: a non-local geometry in which any pair 
of qubits may interact and a local geometry in which the 
qubits are arranged in a one dimensional chain, where 
only nearest neighbor interactions arc allowed. More 
specifically, we used a combination of random single qubit 



rotations and a fixed two-qubit operation, namely the 
controlled phase gate. We found that in this scheme, 
entanglement is generated more quickly in the non-local 
geometry than in the one-dimensional chain. Thus, non- 
local implementations of quantum computation are ex- 
pected to be more efficient in generating highly entangled 
states. 
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